High Fidelity Quantum Gates in the Presence of Dispersion 
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We numerically demonstrate the control of motional degrees of freedom of an ensemble of neutral 
atoms in an optical lattice with a shallow trapping potential. Taking into account the range of 
quasimomenta across different Brillouin zones results in an ensemble whose members effectively 
have inhomogeneous control fields as well as spectrally distinct control Hamiltonians. We present 
an ensemble-averaged optimal control technique that yields high fidelity control pulses, irrespective 
of quasimomentum, with average fidelities above 98%. The resulting controls show a broadband 
spectrum with gate times in the order of several free oscillations to optimize gates with up to 13.2% 
dispersion in the energies from the band structure. This can be seen as a model system for the 
prospects of robust quantum control. This result explores the limits of discretizing a continuous 
ensemble for control theory. 
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I. INTRODUCTION 

Optical lattices [l| are a suitable platform for the sim- 
ulation of strongly correlated many-body and other con- 
densed matter systems 0-11]. These simulations can in 
turn be used to compute otherwise intractable problems 
such as antiferromagnetic transitions and conduc- 
tor to insulator transitions @, 0] j to name a few. Op- 
tical lattices are also a promising candidate system for 
quantum information processing 0. In the simplest 
case, a one-dimensional optical lattice is created by two 
counter-propagating beams of coherent light which cre- 
ate a standing wave. This generates a trapping poten- 
tial in which sufficiently cold atoms can be bound and 
trapped in a one-dimensional periodic potential but can 
move freely in the other two spatial dimensions. It also 
provides an opportunity to investigate a difficult control 
problem in state evolution, particularly from the perspec- 
tive of quantum computation with qubits formed by the 
vibrational states of atoms in the lattice, so-called exter- 
nal qubits [i-fl^. 

One of the strengths of optical lattices is that there is a 
great deal of freedom for control and manipulation, how- 
ever, in practice finding the control sequences to perform 
a desired task can be difficult. The control task we ex- 
amine in this document is to perform a single-qubit gate 
between states in the lowest two Bloch bands. The eigen- 
states of a periodic potential have a band structure, with 
states within the bands distinguished by a quasimomen- 
tum, see Fig. [TJ Our target gates are naturally vertical. 



i.e., they connect states with equal quasimomenta. The 
controls, phase and intensity modulations of the control 
laser, both do not affect the discrete displacement sym- 
metry hence they do not change the quantum number k. 
This can be shown by recognizing that the quasimomenta 
k describe the eigenvalues e**^" of the discrete displace- 
ment operator. The controls we have at our disposal are 
collective controls of the lattice. The challenge is to con- 
trol an ensemble of particles in parallel even though they 
each undergo a different quantum evolution, which math- 
ematically is identical to the problem of robust control 

MM. 



For complex problems such as controlling an ensemble 
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FIG. 1: (color online) The first four energy bands (highlighted 
in grey) for a 1-D optical lattice are shown in the localized ba- 
sis for potential depth r = 13 and compared to the sinusoidal 
trapping potential. Refer to Fig. [2] 
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it is generally difficult to find control sequences analyt- 
ically. An alternative is to use techniques from optimal 
control theory, where the fields are found numerically by 
optimizing a cost function, usually a state or gate fidelity 
j2(]| - |23| . These numerical approaches are very powerful 
because the optimization algorithms they use are mostly 
platform independent and easily extendable to account 
for constraints such as those of ensemble control. Op- 
timal control theory has been applied to the problem of 
robust or ensemble control in many different contexts like 
NMR [23-127}. many-body entanglement [28j . spin-chains 
[29| , and spin systems [3G |. 

The remainder of this paper is a follows. In section HIl 
we introduce the one-dimensional lattice model, its band 
structure and explore the interactions we utilize to con- 
trol this system. In section Hill we develop how this inho- 
mogenous system can be treated within optimal control 
theory. The results of our numerical control treatment 
are shown in section ITVl 



where n indicates the energy band and k the quasimo- 



II. PHYSICS OF THE OPTICAL LATTICE 



A. Model 



mentum and ip^"^ {x + rm:) = ijj^' (x) for all integers m. 
These Blocli functions satisfy the Schrodinger equation 



with the now fc-dependcnt Hamiltonian 



(k) 



cos 2x). 



(2.4) 



(2.5) 



The resulting wavefunctions can be calculated using the 
central matrix method or the Bloch method (33| . Of pri- 
mary importance is that the quasimomentum is a good 
quantum number. The single-particle energy eigenstates 
of the optical lattice thus, from a control perspective, 
form an inhomogenous ensemble of discrete quantum sys- 
tems with different parameters k in the Hamiltonian of 
Eq. (|2.5I) . As we shall see later, the controls also conserve 
fc, thus preserving this ensemble decomposition. 

To control this system we introduce the laser parame- 
ters •q{t) and (j}{t) which represent the ratio of intensity 
with respect to r and the phase, respectively, of the lasers. 
Rewriting (j2.5p in terms of these new variables we obtain. 



We consider a physical model of the optical lat- 
tice similar to the experimental apparatus described in 
Refs [ni HH, mi ■ The optical lattice potential is gener- 
ated by two vertical counter-propagating lasers with an 
incidence angle of 49.6° and is loaded with ^^Rb atoms 
from an optical molasses at 10 /xK. The Hamiltonian for 
an alkali atom interacting with an optical lattice poten- 
tial in 1-D is given by 



Ho = ^ + Uo svc^ikLx) 
Zm 



(2.1) 



where m is the mass of the atom, k^ wave number of the 
lattice laser light and Uq is the lattice potential depth. 
In experiments, the depth is typically 18-30£'fl though 
we consider a wider range for this numerical study. Here 
En = fi?k\/2m is the recoil energy, i.e., the kinetic en- 
ergy the atom gains by absorbing a lattice photon. At 
these lattice depths for typical atomic densities the atoms 
are essentially non-interacting [ll| and thus we can sim- 
plify the model by considering the evolution of a single 
atom. 

The dimensionless form of this Hamiltonian is given 
by, 



Hq ~ P + r sin x ^ p 



-(l-cos2a;), (2.2) 



where r = Ua/Efj is the lattice potential depth in units 
of recoil energies and p is now the momentum in units of 
hk]^. Bloch's theorem states that a periodic Hamiltonian 



such as Eq. 



has eigenstates of the form. 



(fe) 



(2.3) 



H('^\t) ^{p - fc)2 + : [1 + r,{t)] [1 - cos {2x + <j>m 



-H, 



[1 - cos(2x)] 



+ 2 + [1-cos i2x + q^m 
=H^o''^ + ^ [1 - (1 + V{t)) cos (/)(t)] cos(2x) 
+ ^(l + #))sin</)(i) sin(2x) + ^7^0- 

(2.6) 

In order to express the Hamiltonian in the standard form 
for quantum control, (23l |. we rcparameterizc the control 
fields in terms of 



ait) 



T 

m =4 



1 - [1 + 77(t)]cos(/)(t) 
[1 + 77(0] sin(/.(t) 



(2.7) 



so that the total control Hamiltonian (neglecting the 
global phase) is 

H'-'^^t) = H^J"^ +2a{t)cos(2x) + 2p(t)sin{2x). (2.8) 

The lattice will be filled with multiple atoms occupying a 
range of fc- values. Measurement of the success of a quan- 
tum operation can be performed in the manner suggested 
in (b^ I , which involves averaging single particle measure- 
ments over the entire ensemble. Demanding that we find 
a{t) and /3{t) that perform the desired gate irrespective of 
k ensures that any gates we find will retain their fidelity 
after this averaging. 
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B. Dispersion 



C. Relation to the Charge Qubit 



States with different quasimomenta evolve indepen- 
dently under the applied Hamiltonian Hq''^ and so that 
they can be viewed as an ensemble of non-interacting 
systems, labelled by the quasimomentum. The eigen- 
energies of this system depend on the quasimomentum, 



H, 



leading to dispersion of 

the resonance frequencies across different quasimomenta. 
That is the resonant frequencies of excitations associated 
with different vertical transitions are not the same, see 
Fig. [2] We characterize the dispersion by the dimension- 
less quantity 



V 



AE, 



(1) 



AE, 



(0)' 



(2.9) 



'Eq'^'^ . In addition to energy disper- 



where AE^^l^ = e['''' - 
sion there is dispersion in the strengths of the couplings 
between the bands with respect to the quasimomentum. 

The main consequence of the dispersion is that a sin- 
gle resonant harmonic excitation cannot simultaneously 
resonantly excite all transitions. This poses a formidable 
challenge for controlling an ensemble with high disper- 
sion. For the optical lattice, the amount of energy disper- 
sion decreases with an increase in the strength of poten- 
tial. Fig. [21 Unfortunately, the anharmonicity of the lat- 
tice, which is essential to resolve different transition fre- 
quencies, also decreases with increasing potential depth. 
This leads to a trade-off when choosing the depth of the 
lattice between dispersion and anharmonicity. Note that 
the latter alone can be addressed with techniques like 
DRAG 111. 
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FIG. 2: (color online) The first four energy 1-D band struc- 
tures are shown in the momentum basis for four different po- 
tential depths r. The lower potential depth r = 2 (top left) 
shows the largest amount of energy dispersion with the clos- 
est energy crossings. The depth r = 30 (bottom right) shows 
a lesser amount of dispersion and larger energy splittings. 



A device that behaves similarly to the optical lattice 
in the context of superconducting qubits [ssl l36j is the 
charge qubit [37l - l40j . The effective Hamiltonian of this 
system takes the form 



4Ecin - rig)^ - Ej cost 



(2.10) 



where n is the number of Cooper pairs on the supercon- 
ducting islands, (j) is the phase difference to the island, 
Ec and Ej are the charging and Josephson energy, re- 
spectively and Ug is the effective offset charge. 

Since n and (f> are canonically conjugate variables, we 
can directly map this Hamiltonian onto that of the opti- 
cal lattice by letting 



Ej r 

n ^ p, (j) ^ X and — > -. (2-11) 

Ec 2 

ate charge, 
. The gate 

charge takes on a single value, unlike the the quasimo- 
mentum, but has uncertainty due to noise from fluctu- 
ating charges and voltages These fluctuations are 
slow and in particular they occur on time scales com- 
parable to the duration of an experiment. An ensem- 
ble in these experiments, which are carried out on single 
quantum systems, is built up by repeating the exper- 
iment. In order to control the charge qubit we must 
account for this inherent spread over Ug. The solution 
to this problem for the transmon qubit is to work in 
a lower dispersion reg ime, with a tradeoff in the effec- 
tive anharmonicitv|39l |40| . The optimization techniques 
for these systems will be identical to the optical lattice, 
even though the causes of the dispersion are very differ- 
ent: ensemble averaging versus parameter uncertainty. 
Note that in superconducting qubits, the CORPSE pulse 
sequence [42| has been applied to mitigate fluctuations of 
the energy splittings [43l|. which has been investigated 
further numerically |44| . 



The only remaining component is to let the 
rig, play the role of the quasimomentum 



III. TECHNIQUE FOR OPTIMIZATION 

In order to find pulses that prepare a desired state or 
unitary transfer we use numerical techniques from opti- 
mal control [lO, [2^ I45l - l47j . Primarily we use a variant 
of the GRAPE (gradient-ascent pulse engineering) algo- 
rithm to numerically find optimal control fields for a fixed 
duration pulse [2^ . GRAPE is a numerical optimization 
method that maximizes the fidelity, $, with respect to 
the control fields, u, for generating gates or states from 
some initial configuration of the system. For the major- 
ity of this section we will focus on gate preparation since 
it is of greater interest and because the state preparation 
protocol is a straightforward modification. 

The optimization follows the form of a gradient search, 
where at each iteration, j, the control fields (initialized 
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to a value u^^^ ) arc updated according to 

u(j'+i) = u(^) + eV„$, (3.1) 

where e is a small step parameter, usually chosen adap- 
tively for faster convergence. An essential feature of an 
optimal control algorithm such as GRAPE is to compute 
the gradient V„$ in an efficient manner [2^|4^,|43. This 
protocol is not guaranteed to find a global optimum; how- 
ever for many problems it is easy to find high-fidelity so- 
lutions after optimizing a small number of initial guesses 

IMIiii- 

We require a suitable fidelity function to average the 
Hamiltonian over the ensemble. As a starting point we 
take the standard average gate fidelity, 



^|Tr(ft[/)|2 



(3.2) 



where V is the desired gate that we wish to implement, 
U is the total time-evolution due to our choice of control 
fields and d is the dimension of the Hilbert space. Our 



desired gate should act on the 



states independent 



dkVIl{k), where n(fc) is 



of fc, so we can write it as V 
the projector onto the quasimomenta in the first Brillouin 
zone and V is an operator that acts only on the energy 
band degree of freedom, thus commuting with JI{k). If 
we define 



(3.3) 



then the total evolution operator of our system is U 
/dfc'C7('=')n(fc'), leading to the fidelity 



d2 
1 

1p 



Tr (J I dkdk' V'^Tl{k)U^^"^Tl{k' 



(3.4) 



where the final trace is taken over only the band degree 
of freedom. 

This fidelity is not the same as simply averaging the 
standard gate fidelity over the ensemble. An important 
feature of Eq. p.4|) is that the absolute value is taken 
after the ensemble average, forcing the phases for differ- 
ent members of the ensemble to be equivalent in order 
to maximize the fidelity. This global phase on each U'^^^ 
at first seems like it should be irrelevant, but it plays an 
important role when it comes to optimizing pulses with 
a finite bandwidth. A simple example where this global 
phase is important is the following: restrict to evolutions 
in SU{d) where all Hamiltonians are traceless and all 
unitary evolutions therefore have det(?7) = 1. In this 
setting there is a finite set of unitary maps that maxi- 
mize I Tr(Ftc/)| namely U = e''^'"i''^V, where q = l...d 
and we assume dei{V) ~ 1. These are the only global 
phases for which U has a unit determinant. Ignoring this 
phase and attempting to optimize a gate that is invariant 



to a parameter such as the quasimomentum can lead to 
a situation where the global phase locks to different val- 
ues in different regions of fc-space. The finite bandwidth 
of our controls does not allow U'^^^ to change instanta- 
neously with k, and so there will be some intermediate 
regimes between two values of the global phase where the 
gate performs very poorly. The evolution of the optical 
lattice we have outlined in this paper (Eq. 12. 8p is not 
restricted to SU{d), though we suspect that a similar ar- 
gument can be made. This requires further investigation. 
In practice our numerical optimization can get trapped if 
we neglect this global phase. We see that there are high 
fidelity regions in k that are connected by intermediate 
regions with low fidelity. Our fidelity accounts for this 
issue by insisting the global phase to be uniform across 
the ensemble. 

In order to make this fidelity tractable for numerical 
evaluation we sample the quasimomentum for a discrete 
number of values ki,k2, ■ ■ ■ kM and then we arrive at the 
fidelity 



{dAiy 



M 



^Trfytf/C^^ 



(3.5) 



The Nyquist sampling rate sets a minimum bandwidth 
of the control pulses to ensure a minimum fidelity with 
respect to variations in the quasimomenta. Therefore if 
we sample the quasimomentum on a fine enough grid we 
can be assured that the intermediate values have a high 
fidelity, and as verification we can test this hypothesis 
after optimizing our fields. 

For GRAPE we must provide the derivative, d^/duq, 
which we calculate by 



9$ 



duq {dMy 



M 



-Re^Tr 



duo 



M 



(3.6) 



i'=i 



Essentially, we must calculate 817^'''^ duq for M different 
evolutions leading to a computational-scaling for evaluat- 
ing this derivative that is only M times more than stan- 
dard GRAPE on a d-dimensional system. However, this 
is significantly faster than optimizing a gate for a dM- 
dimensional Hilbert space since multiplying dxd matrices 
scales as 0{d^). 



IV. RESULTS 

Using the techniques outlined in the previous section, 
we optimize control parameters for the optical lattice 
Hamiltonian in Eq. (|2.2p for varying potential depths 
and pulse durations. Our target gate is an X^r-gate on 
the first two energy bands. We consider both a Rabi 
pulse assuming no dispersion and bounded random con- 
trols fields as the initial guesses for the control fields. 
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FIG. 4: (color online) Comparison of gate error for well- 
behaved optimized controls in Fig. [3] where the blue dotted 
line corresponds to control sequence (a) and the red solid line 
to (b). For (a) the gate error is at most 1.8% across the pos- 
sible quasimomenta and the average gate error is 0.6%. In 
(b) the gate error has been minimized for the quasimomenta 
points that were sampled for optimization but was 66.8% on 
average when sampled over 100 points. 
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FIG. 3: (color online) Comparison of two optimized con- 
trols for preparing a X^-gate when the optimization uses only 
ten points in quasi momentum space. Both pulses are well- 
behaved in the sense that they don't translate the lattice by 
a full lattice site or possess very large excursions from the ini- 
tial value for r. (a) An optimized pulse for a potential depth 
r = 17 with 5.4% dispersion and fidelity 99.3% which is cal- 
culated by sampling over 100 quasi momentum values. The 
duration of the pulse was 5 free oscillations (at k=0). (b) 
An optimized pulse for potential depth r = 12 and disper- 
sion 13.2% with a sampled fidelity of only 66.8% . The time 
duration of the pulse was 25 free oscillations (at k=0). 



The evolution is simulated using only the first six energy 
bands from our model and sampling over 20 values of 
the quasimomenta. While we sample over only a small 
number of values for the quasimomenta in order to have 
efficiently performed the numerical search, when we cal- 
culated the final fidelity we sampled over a finer grid in 
quasimomentum space. 

Here are two examples of optimized controls for per- 
forming ATTT-gates on our optical lattice model. One is a 
control for an optical lattice system with potential depth 



r = 17 (which corresponds to a dispersion of 5.4%), as 
shown in Fig. 3(a) The duration of this control was 5 
free oscillations (at k=0). The gate error across the en- 
semble was less than 1% with the exception of particles 
with quasimomentum near zero. Fig. 21 Even with such 
high dispersions in the energies and control Hamiltonians 
we are able to find gates with reasonable fidelities. 

The second control in Fig. |3(b)| is for an optical lattice 
system with potential depth r = 12 and thus a dispersion 
of 13.2%. The average fidelity for the optimized points 
was 99.96% but after finer examination of the fidelity 
across the quasimomenta it was found to have an aver- 
age of 66.8%, as shown in Fig. Thus a more cautious 
approach should be taken with regards to optimizing a 
coarse sampling over the quasimomenta space. These 
limitations can be overcome by finer sampling over the 
range of quasimomenta, however, this can become com- 
putationally expensive. 

We observe the relationship between the maximum fi- 
delity of the solutions and the level of dispersion of the 
Hamiltonian. Ideally, our gradient search will halt when 
a local maximum in the fidelity is reached but this may 
take a considerable amount of computational resources. 
All our optimizations were halted after either the algo- 
rithm converged to a solution or 10^ updates in the con- 
trol parameters were performed. Nevertheless, for short 
times (Fig. [S]) we find excellent control fields ($ > 0.90) 
for a range of potential depths from 12 < r < 110. 

In general, for long gate times, the maximum fidelity 
solution for X^-gates became lower as the dispersion 
was increased. Fig. [51 The optimization becomes less 
tractable with higher dispersion since we are asking a 
single set of control fields to solve substantially different 
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problems depending on where the system is in quasimo- 
mentum space. As a result, broadband pulses must be 
tailored to accommodate a range of possible energies and 
couplings. 

We also find that for high dispersion the fidelity be- 
comes worse as the gate time becomes longer. This may 
seem counterintuitive, but it is due to the fact that as 
the gate time increases the fidelity can vary more quickly 
as a function of k. For very short gate times, the fidelity 
is fairly constant across quasimomentum space, but for 
long gates we find that the fidelity is only high for the 
specific points we optimized, and dips to almost zero in 
the intermediate regimes, as seen in Fig. |4l 

We can see the connection between gate duration and 
rate of change of fidelity very easily through the simula- 
tion described in Fig. [T] Here we run GRAPE for exactly 
one value of the quasimomentum and and then look at 
the performance of this control field across k. We find 
that as the duration of the gate increases the fidelity func- 
tion becomes more tightly peaked about the single value 
we have optimized. In principle, there should alway be 
higher fidelity control fields at longer gate times but find- 
ing these fields becomes computationally expensive as we 
need to sample more values in quasimomentum space or 
simply run the algorithm from more initial conditions in 
the hope that one produces a flat fidelity curve. 
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FIG. 6: (color online) Same as Fig. [5] but witli times from 15 
to 70 free oscillations (at k=0). 
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FIG. 7: (color online) Pulses for performing X^-gate at r = 2 
for three different times that were optimized specifically for 
k=0.5 and the fidelity response over quasimomenta shown. 
This shows the effect of Nyquist bandwidth limits of the pulses 
on fidelity. Shorter pulses have larger spectral bandwidth and 
thus affect a larger range of quasimomenta. 



FIG. 5: (color online) The maximum fidelity for an optimized 
pulse over a range of times, from 1 to 10 free oscillations (at 
k=0), preparing a X.^-gate as a series of different potential 
depths. The potential depth ranged from r = 0.25 to r = 110, 
giving a dispersion ranging from 93.8% to 0.01%, respectively. 
Each point is an average of optimized pulses from 1 Rabi and 
10 random initial pulses. 



V. CONCLUSION 

We simulated our controls to observe the fidelity for 
all possible quasimomenta in the ensemble. We also ob- 
served the minimum time to prepare a A"7r-gate for a 



given dispersion using our techniques. This shows that it 
is possible to perform operations of high fidelity in rea- 
sonable times for a system with inhomogeneity in both 
energies and couplings. Specifically, we have demon- 
strated control of a shallow one-dimensional optical lat- 
tice with atoms that have inhomogeneous control matrix 
elements and distinct spectra. With an ensemble aver- 
aged optimal control technique we have demonstrated ro- 
bustness and shown that with 13.2% dispersion in ener- 
gies it is possible to achieve fidelities as high as 98.3% 
with relatively short pulses. In lower potential depths, 
some of the added dispersion can be corrected allowing 
to take advantage of greater anharmonicity to perform 
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shorter gates. 

The technique does contain Umitations due to coarse 
sampling of the ensemble for optimization. With finer im- 
provements, this should be resolved. We have also shown 
that this technique performs better for shorter times with 
high dispersion due to the flatter response of the fidelity 
across the ensemble. 

Where developing control pulses from analytic tech- 
niques can become intractable for complicated systems, 
optimal control provides a way out. The control tech- 
niques outlined in this paper are not exclusive to the 
optical lattice and can be applied to other quantum sys- 
tems. In fact, these techniques can be applied to any 
system with noise in control parameters, energy fluctu- 
ations, ensembles with varying couplings and energies, 
or of uncertain system specifications. With further en- 
hancement of numerical and analytical techniques, there 
is no reason that higher fidelities and shorter times should 
be impossible. Much work has been done in this field 
however there is potential for improved techniques for 
designing robust control pulses, allowing for greater con- 
trol over quantum systems [soj . What has been shown in 



this paper is the application of optimal control to robust 
pulse design. 
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